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ABSTRACT 



Aims. The supergiant fast X-ray transient source IGR J 16479-45 14 was observed in outburst two times with Swift. Its quiescent state 
was investigated in-depth only once in 2008 through a relatively long pointed observation with XMM-Newton. The latter observation 
was taken about 1.7 days after the outburst in 2008, and showed an X-ray eclipse-like event, likely caused by the supergiant compan- 
ion. At present, this is the only supergiant fast X-ray transient that displayed an evidence for an X-ray eclipse. 

Metliods. Here we carry out a comparison between the most recent outburst of IGR J 16479-45 14, caught by Swift on 29 January 2009 
and those detected previously from this source. 

Results. The decay from the outbursts in 2005, 2008 and 2009 presents many similarities, and suggests a common mechanism that 
modulates the mass accretion rate onto the neutron star in IGR J16479-4514. 

Key words. X-rays: binaries - binaries: eclipsing - stars: individual (IGRJ16479-4514) -stars: neutron - X-rays: stars 



1. Introduction 

IGR J16479-4514 is one of the supergiant fast X-ray transients 
(SFXT), a subclass of high mass X-ray binaries (HMXBs) host- 
ing an accreting neutron star (NS) and a supergiant companion 
(Sguera et al., 2006; Walter et al., 2006). These sources undergo 



and thus very little is known on the mechanism that triggers the 
transitions between outburst and quiescence. 

IGR J 16479-45 14 was observed w ith Swift in outburst two 
tim es, on 2005 August 30^guera et a l., 2008) and 2008 March 



sporadic outbursts (peak luminosities of -10^^''- 10"" erg s"'), 
with typical durations of few thousands of seconds, and remain 
in quiescence for most of the time, with a mean X-ray luminosity 
in the -10^^-10^''* erg s"' range. In order to explain the behaviour 
of the SFXT sources, a number of different scenarios have 



been proposed din'tZandl 120051: loiter & Zuritall eras, 2007; 
iNegueruela etall l2008t ISidoh et all 120071: iBozzo et al.. ,2008i) . 



All these models involve a NS accreting from the intense wind 
of its supergiant companion, and interpret the transitions from 
outburst to quiescence as being due to strong variations in the 
mass accretion rate onto the NS. These variations have been as- 
cribed to a centrifug al and/or magnetic gating mechanism by the 
NS magnetosphere dBozzo et al., 2008]), or to clumps and struc- 
tures in the wind of the supergiant star ([W alter & Zurit a Herai, 
120071; iNegueruela et al.Ll2008h . 1sidoli et aLi (2007.) proposed that 
SFXT outbursts result from a neutron star in a tilted orbit cross- 
ing a dense equatorial disk around the supergiant companion. As 
different models make different predictions with respect to the 
origin of the quiescent emission in SFXTs, an in-depth knowl- 
edge of the timing and spectral properties of the low activity 
states of these sources is expected to help discriminating be- 
tween the above models (iBozzo et al.L 12008). At present, only 
a few observations of quiescent SFXTs have been carried out, 



19 ( Romano et al.L l2008 a'!?). Its X-rav emission properties dur 
ing these bright events is fairly well known. On the contrary, the 
quiescent state of IGR J16479-4514 was investigated in-depth 
only once with XMM-Newton^ This observation lasted ~32 ks 
and was carried out about 1.7 days after the bright outburst in 
2008, with the aim of gaining insight on the mechanism that 
drive the quiescent emission of this source. Surprisingly, XMM- 
Newton caught IGR J 16479-45 14 undergoing a transition from a 
bright state (2-10 keV absorbed flux 1.0x10"" erg/cm^/s) to a 
lower flux state (2-10 kev absorbed flux 7.5x10"'-' erg/cm^/s). 
The spectral and timing analysis of this observation showed un- 
ambiguously that the transition was due to the obscuration by 
an occulting body, most likely an eclipse by the supergiant com- 
panion. This revealed that at least some of the very large X-ray 
luminosity swings displayed by the SFXT sources might not be 
due only to genuine changes in the mass accretion rate. 

Here we report on the most recent outburst of IGR J 16479- 
4514, detected by Swift on 2009 January 29, and compare this 
outburst with those observed previously in 2005, 2008. 



2. Data analysis and results 

Swift /BAT caught a new bright o utburst from IGR J16479-4514 
on 2009 January 29 dRomano et al.. 2009). The start time of the 
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' IWalter et aP d2006) also reported on a -2.9 ks XMM-Newton obser- 
vation of this source in quiescence, but it was alfected by a very strong 
and variable background. 
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Fig. 1. BAT and XRT/PC light curve of the 2009 outburst (time 
bin is 50 s). 

BAT trigger was 06:33:07 UTC (including UTCF correction). 
Swift /XRT slewed to the source about 800 s after the BAT trig- 
ger, and observed in Windowed Timing (WT) mode during the 
first 46 s. Data in photon counting (PC) mode were then accu- 
mulated for a total exposure time of 2.6 ks. We analyzed the data 
from observation ID. 00341452000, that includes the BAT data 
of the outburst, as well as the following XRT WT and PC data, 
and use also Swift follow-up observations that were carried out 
up to 20 days after the BAT trigger. 

In order to compare the 2009 outburst with those detected 
previously, we also analyzed Swift observations of this source 
that were carried out during the outbursts in 2005 an d 2008 
(these observations w ere previously published by Sguer a et al.L 
120081: iRomano et"an l2008a) . We used outburst and follow-up 
observations up to ~20 days after each BAT trigger. A log of the 
data sets is given in Table [1] 

The Swift data were analyzed by using standard procedures 
(iBurrows et al .. 2005) and the latest calibration files available. 
The BAT and XRT data were processed with the batproduct 
(v.2.42) and xrtpipeline (v.0.12.1) tasks, respectively. Filtering 
and screening criteria were applied by using ftools (Heasoft 
v.6.6.1). We extracted source and background light curves and 
spectra by selecting event grades of 0-2 and 0-12, respectively 
for the WT and PC mode. Exposure maps were created through 
the XRTEXPOMAP task, and we used the latest spectral redis- 
tribution matrices in the heasarc calibration database (v.Oll). 
Ancillary response files, accounting for different extraction re- 
gions, vignetting and PSF corrections, were generated by using 
the xrtmkarf task. When required, we corrected PC data for 
pile-up, and used the xrtlccorr to account for this correction 
in the background-subtracted light curves. 

For each observation in Table [1] we extracted the light curve 
and spectrum, and derived a mean X-ray flux by fitting this 
spectrum with an absorbed power law model (we used Xspec 
V. 12.5.0). Spectra with an adequate statistics were rebinned in or- 
der to have at least 20 photons per bin so as to permit fitting. 
Lower statistic spectra were rebinned in order to have a t least 
5 pho tons per bin and then fit by using the C-statistics (ICashL 
Il979h . Observations with less than 100 counts were treated sep- 
arately, as the C-statistics is not recommended with such a small 
number of photonfl In these cases we estimated the source 
count rate of the observation with sosta (ximage V.4.4.1), and 

^ https://astrophysics.gsfc.nasa.gov/XSPECwiki/statisticaLmethods 
JnJfSPEC. 
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Fig. 2. BAT and XRT/PC spectra of the 2009 outburst. The BAT 
spectrum is integrated over the total outburst duration, whereas 
the XRT/PC spectrum is integrated over the A-i-B+C time inter- 
vals shown in Fig. [T] The former spectrum was fit by using a 
simple power law, whereas for the latter we also included the ef- 
fect of the absorption at low energies (the best fit parameters are 
reported in Table [T]!. The lower panel shows the residuals from 
the best fits. 

then used this count rate within webpimm^ in order to derive the 
X-ray flux (we assumed the same spectral model of the closest 
observation for which a spectral analysis could be carried out). 

Figure [U shows the BAT (15-25 keV) and XRT/PC (0.3- 
10 keV) light curve of the 2009 outburst (XRT/WT data in the 
0.3-10 keV energy band are included only in Fig. As the 
observation interval of BAT and XRT/PC do not significantly 
overlap, we did not attempt to fit a combined BAT+XRT spec- 
trum. Instead we fit them separately and report the results of 
these fits in Table [1] The BAT spectrum, integrated over the 
total duration of the outburst, and the XRT/PC spectrum of 
observation ID 00341452000 (A+B-t-C intervals) ai-e shown in 
Fig. |2] together with the best fits and the residuals from these 
fits. We also searched for spectral variability by extracting time 
resolved spectra in the intervals A, B, and C of Fig. [1] This 
analysis revealed that the source X-ray flux decreased from 
3.2x10"'° erg/cm^s"' in interval A, to 3.6x10"" erg/cm^/s"' in 
interval C. The best fit powerlaw index and the absorption col- 
umn density did not show any significant variation (they were all 
constant within the errors). 

In Fig. [3] we plot the evolution (up to -20 days from the 
BAT trigger) of the absorbed 2-10 keV X-ray flux during the 
2005, 2008 and 2009 outburst of IGR J 16479-45 14. In this figure 
we used the fluxes from Table [1] and considered half of each 
observation duration as the error on the time axis. 

Finally, in Fig. |4] we show the combined XRT/PC and 
XRT/WT light curves of the three outburts in detail. For clarity, 
we divided the count rate of the first 5 points (XRT/WT data) of 
the 2008 light curve by a factor of 10. 

It is seen from Fig. [T]that the rise and decay times of the 
2009 outburst are of order of a few thousand seconds, simi- 
lar to the time scales measured in the 200 5 and 2008 outbursts 
(ISguera et all 120081: iRomano et al.Ll2008ah . The fluxes and spec- 
tral properties of the source emission during the peak and the 
beginning of the outburst decay (up to ~7 ks after the BAT trig- 
ger) are also in agreement with previously reported results. The 
typical outburst spectrum at high energy (15-50 keV) can be fit 

^ http://heasarc.nasa.gov/Tools/w3pimms.html 
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Table 1. Observations log of IGRJ16479-4514 (see note at the bottom of the Table for other specifications). 
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NOTE: Spectra extracted from these observations are fit with an absorbed power law (absorption column density Nh and photon index F). Fobs is the 
XRT/PC or XRT/WT (BAT) absorbed flux in the 0.3-10 keV (15-80 keV) energy band. EXP indicates die total exposure time of each obsei-vation 
( Swift observations comprise several snapshots and are not continuous pointings at the source). 
" 3o" upper limit. 



' C-statistic has been used to fit the spectrum. 
In this case the BAT flux is in the energy range 15-50 keV. 



by a power law of photon inde x ~2.5. The mean X -ray flux is 
of few times 10"'' erg/cm^/s"' (ISguera et al.L |2008|) . The 2-10 
keV spectrum is flatter and can be well fit by using a power 
law of photon i ndex r~l and an abso rption column density of 
-5x10^^ cm-^ (iRomano et all 12008 a). We did not attempt to fit 
the broad-band spectrum of the 2009 outburst, since the obser- 
vation intervals of BAT and XRT do not overlap significantly 
(see Fig. Moreover, owing to poor statistics of the XRT/PC 
data, time resolved spectra could not be investigated during the 
first 7 ks after the BAT trigger (intervals A, B and C in Fig. [TJ. 
Therefore, we could not search for spectral variations during the 
decay from the outburst. However, we note that follow up obser- 
vations carried out up to 20 days after the BAT trigger showed a 
moderate increase of the absorption column density (a factor of 



~2) and a steepening of the power law index (from ~ 1 to ~2, see 
Table[B. 

In order to compare the 2009 outburst with those detected 
previously in 2005 and 2008, we studied the variation of the X- 
ray flux from IGR J 16479-45 14 during the 20 days that followed 
each of these outbursts. This is shown in Fig. [3] In the middle 
panel of this figure we give the X-ray flux of the 2008 outburst; 
the horizontal dotted lines represent the high and lower flux 
state of IGRJ16479-4514 as observed in 2008 hy X MM-Newton 
see Se ct. [T]l- According to the results discussed in lBozzo et al.l 
2008bh . these correspond to the average X-ray flux during the 
eclipse ingress and the mid-eclipse, respectively. For compar- 
ison, these two dotted lines are also plotted in the upper and 
lower panels of Fig. [3] that give the evolution of the X-ray flux 



4 



E. Bozzo et al.: IGR J 16479-45 14 in outburst (RN) 



□ 


ffi. . 




-T- 


2009 : 

-1 1 








# 









5 


10 


15 


30 


' a 

r ffi 


f 






200B 1 









5 


10 


15 


30 


□ 








2005 - 









5 


10 


15 


30 



days after BAT trigger 

Fig. 3. Evolution of the X-ray flux during the 2005, 2008, and 
2009 outbursts. We used all the observations reported in Table[Tl 
and considered the BAT trigger time as the origin of the X axis. 
The BAT ti'igger times ai-e 2005-08-30 04:08:48 UTC, 2008-03- 
19 22:44:45 UTC, and 2009-01-29 06:33:07 UTC. In all panels, 
squares indicates the BAT X-ray flux (15-80 keV for the 2009 
and 2008 outburst and 15-50 keV for the 2005 outburst; we did 
not indicated the errors for these fluxes), open circles and trian- 
gles are for XRT/PC and XRT/WT fluxes (absorbed flux in the 
2-10 keV band), respectively. The arrows indicate 3cr upper lim- 
its. Errors on the X axis correspond to the half duration of each 
observation. 



during the 2009 and 2005 outbursts, respectively. We note that, 
in all the three cases, several observations have an X-ray flux 
that falls inside the range defined by the two dotted lines. As 
a consequence of this and the relatively sparse time coverage, 
the plots in Fig. [3] do not permit to single out the presence of 
an X-ray eclipse. Rather IGR J 1 6479-45 1 4 may simply alternate 
between higher and lower flux states; unfortunately, the origin of 
these states can hardly be investigated without a detailed spectral 
ana lysifl 

By setting the start times of all light curves in Fig. |3]at the 
onset of the outbursts, we can compare the flux evolution and, 
in particular, verify whether the eclipse-like event ~1.7 d after 
the 2008 event, has a counterpart in the 2005 and 2009 light 
curves for comparable delays. Stated differently, this amounts 
to checking whether the outburst and the eclipse-like event are 
driven by the same, possibly orbital, "clock". 

Swift observations in 2008 clearly show that the source X-ray 
flux between 1-2 days after the BAT trigger was compatible with 
the lower flux measured by XMM-Newton (i.e. the mid-eclipse 
state). Furthermore, the mean X-ray flux of the subsequent ob- 
servation in 2008 (that is compatible with the eclipse egress, if 
the eclipse is symmetric), suggests an upper limit to the duration 
the eclipse-like event of about ~3 days. Even thought the orbital 
period of IGR J16479-4514 is still not known, this limit seems 

The 2008 XMM-Newton spectra, on the contrary, provided convinc- 
ing evidence that the lower intensity state resulted from an occultation 
episode. 



reasonable, as i n other SFXTs ari orbital period of ~30 days ha s 
been measured (iBird et all l2009t: IZurita-Heras & Chatvl l2009l) . 
and, of course, an eclipse can last up to half of the orbital period. 
However, if we supposed that the eclipse-like event takes place 
at a fixed delay from the beginning of the outburst, then the 2009 
observation would indicate that the eclipse duration is shorter. In 
fact, from the upper panel of Fig. [3] we note that the source was 
observed at an X-ray flux comparable with the eclipse egress 
about ~2 days after the BAT trigger, thus limiting the eclipse 
duration to less than one day. A comparison with the 2005 out- 
burst is more complicated because the XRT observation caiTied 
out after the BAT trigger has a large gap close to the expected 
eclipse ingress and egress (see below). 

In order to further study the source behavior during the 2005, 
2008, and 2009 outbursts, we also investigated in detail the X- 
ray light curves of all observations in Table [1] This is shown in 
Fig.|4](times on the X axis are measured from the start time of the 
BAT trigger). In the middle panel of this figure we represented 
with dotted vertical lines the time interval in which the 2008 
eclipse-like event was seen. Under the hypothesis that the X-ray 
eclipse of IGRJ16479-4514always follows the outburst, the two 
dotted lines in the upper and lower panel of Fig.|4]correspond to 
the position where an eclipse should have been observed after 
the outbursts in 2009 and 2005, respectively. Even if the poor 
statistics of the Swift data did not permit a detailed spectral anal- 
ysis, and thus did not provide further support in favour of this 
hypothesis, we note that the source count-rate at these positions 
is consistent within all the three outbursts. 

Apart from the X-ray eclipse-like event, we also note from 
Fig. |3] that the three outbursts displayed many other similar- 
ities. Particularly interesting is the fact the source showed in 
all the reported outbursts an increase in the X-ray activity af- 
ter about 6 ks from the BAT trigger (note that, close to this 
time interval the higher source count rate triggered XRT obser- 
vations in WT mode in all three cases). Furthermore, we note 
the presence of a bright flare (XRT/PC count rate > 1 cts/s) about 
~10 days after the B AT triggers in 2009 and 2005 (see also, 
iLa Parola et aLl l2009l) . The presence of this flare could not be 
investigated during the outburst in 2008, because no X-ray ob- 
servation of IGR J16479-4514 was caiTied out at this period. 

3. Discussion and Conclusion 

The similarities we pointed out in the post-outburst decay of 
IGR J16479-4514 in 2005, 2008 and 2009 suggest the presence 
of a stable structure in the system (at least for a few years), that is 
able to modulate the mass accretion rate onto the NS in a fairly 
reproducible fashion and gives rise to the observed X-ray ac- 
tivity. At present, it is not clear if other SFXT display a similar 
post-outburst behaviour. If this was the case, the implications for 
SFXT models would have to be assessed. 

After this paper was submitted for publication, iJain et all 
(2Q09) reported the discovery of 0.6 d long X-ray eclipse and 
a 3.3194 d orbital period in the Swift/BAT data of IGRJ 16479- 
4514. These result s, besides confirming the interpretation by 
iBozzo et"aLl (l2008h that the sudden drop-off in the source X- 
ray flux detected by XMM-Newton in 2008 was due to an X- 
ray eclipse, provided also a clear confirmation of the conjecture 
proposed in the present research note, namely that the source 
outbursts of 2005, 2008 and 2009 preceded the X-ray eclipse (in 
2008) and the lowest X-ray fluxes intervals (in 2005 and 2009) 
by the same amount of time. Inde e d by u sing the orbital period 
and eclipse duration of iJain et al.l (l2009l) and the eclipse onset 
time from the XMM-Newton observation of 2008 plus 0.3 d to 
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derive the eclipse center time (leading to a superior conjunction 
epoch of MJD 54547,05418, which we adopt as orbital phase 
0), we verified that the onset of the 2005, 2008, and 2009 out- 
bursts, as measured by the Swift /BAT trigger time, occurred at 
an orbital phase of 0.36+0.09, 0.3654+0.0002, and 0.36+0.03 
respective l y. Mo reover the 0.6 d eclipse duration determined by 
Ijain et al] (12009^ agrees well with what we suggested in Sect. |2] 
based on the Swift /XRT data of the 2009 outburst. 

These findings provide strong evidence that the source out- 
bursts, though rare, take place around an orbital phase of 
~0.4 and are thus inherently connected to the orbital motion, 
though the conditions that lead to their occurrence are met only 
rarely. Outbursts recurring at the orbital period were observed 
in three other SFXTs (SAX J18 1 8.6-1703, IGR Jl 1215-5952 
and IGR J 1 8483-031 l:'Bird et a ll l20q9HZurita-Heras & Chatvl 
[2009; Si doli et al.L [2007; Sguera et al.L |2007|) . and are reminis- 
cent of Type I outbursts in Be star X-ray transients. These 
episodes of intense accretion are likely related to periastron 
passage, where the supergiants' wind dens ity is higher and it s 
velocity relative to the neutron star lower dStella et al.L Il986l) . 
Alternatively the outbursts might be due to the neutron star 
crossing an equatorial dense disk ejected by the supergiant 
(ISidoUet all 120071 
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Fig. 4. XRT/PC and XRT/WT light curves of the three outbursts 
in 2005, 2008, and 2009 (time bin is 200 s). All the reported 
times are measured from the BAT triggers (see the caption of 
Fig. |3]l. In the light curve of the 2008 outburst, the beginning 
and the end of the XMM-Newton observation that could have in- 
dicated an X-ray eclipse are marked with dotted vertical lines. In 
the light curves of the other two outbursts, the dotted lines indi- 
cate where the suspected eclipse should have been observed. For 
clarity, we divided the count rate of the first 5 points (XRT/WT 
data) of the 2008 light curve by a factor of 10 (these points are 
separated from the others through a dashed line). 



